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In a generic dynamical system chaos and regular motion coexist side by side, in different 
parts of the phase space. The border between these, where trajectories are neither unstable 
nor stable but of marginal stability, manifests itself through intermittency, dynamics where 
long periods of nearly regular motions are interrupted by irregular chaotic bursts. We dis- 
cuss the Perron-Frobenius operator formalism for such systems, and show by means of a 
1-dimensional intermittent map that intermittency induces branch cuts in dynamical zeta 
functions. Marginality leads to long-time dynamical correlations, in contrast to the expo- 
nentially fast decorrelations of purely chaotic dynamics. We apply the periodic orbit theory 
to quantitative characterization of the associated power-law decays. 

§1. Introduction 

In fluid dynamics one often observes long periods of regular dynamics (lami- 
nar phases) interrupted by irregular chaotic bursts, with the distribution of laminar 
phase intervals well described by a power law. This phenomenon is called inter- 
mittency 1 ^' 2 ^, and it is a very general aspect of dynamics, a shadow cast by non- 
hyperbolic and marginally stable phase space regions. Complete hyperbolicity is 
indeed the exception rather than the rule: almost any dynamical system of interest 
exhibits a mixed phase space where islands of stability or regular regions coexist with 
hyperbolic regions with dynamics mixing exponentially fast. The trajectories on the 
border between chaos and regular dynamics are marginally stable, and trajectories 
from chaotic regions which come close to marginally stable regions can stay 'glued' 
there for arbitrarily long times. These intervals of regular motion are interrupted 
by irregular bursts as the trajectory is re-injected into the chaotic part of the phase 
space. How the trajectories are precisely 'glued' to the marginally stable region is of- 
ten hard to describe. What coarsely looks like a border of an island in a Hamiltonian 
system with a mixed phase space, will under magnification dissolve into infinities of 
island chains of decreasing sizes and cantori 3 ) . 

The existence of marginal or nearly marginal orbits is due to incomplete intersec- 
tions of stable and unstable manifolds in a Smale horseshoe type dynamics. Following 
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the stretching and folding of the invariant manifolds in time one will inevitably find 
phase space points at which the stable and unstable manifolds are almost or exactly 
tangential to each other, implying non-exponential separation of nearby points in 
phase space or, in other words, marginal stability. 

How to deal with the full complexity of a typical Hamiltonian system with 
a mixed phase space is a very difficult and still open problem. Nevertheless, it 
is possible to learn quite a bit about intermittency by considering rather simple 
examples. Here we shall restrict our considerations to 1-dimensional maps which in 
the neighborhood of a single marginally stable fixed point at x=0 take the form 

x f(x) ~ x + cx 1+s (1-1) 

for x — > and are expanding everywhere else. Such a map may allow for escape or 
the dynamics may be bounded, like the Farey map 

T ^ fM -I fo(x) = x/(l-x) x€ Mo = [0,1/2] , , 

I( '~ I fi{x) = {l-x)/x x G M\ = (1/2, 1] ' [L ' Z) 

which has been studied intensively in connection with continued fraction expansions 
and circle maps 4 ) -11 ). 




0.5 





Fig. 1. (a) A tent map trajectory, (b) A Farey map trajectory. 

The dynamics of the Farey map, which is of the form (1-1) with s = 1, is 
compared to that of the hyperbolic, piecewise linear tent map 

fo(x) = 2x xeM = [0, 1/2] 
f 1 (x) = 2(l-x) x G Mi = (1/2, 1] 



x 



fix) 



(1-3) 
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in fig. 1. In sharp contrast to the uniformly chaotic trajectory of the tent map, the 
Farey map trajectory alternates intermittently between slow regular motion close 
to the marginally stable fixed point, and chaotic bursts. The presence of marginal 
stability has striking dynamical consequences: correlation decay may exhibit long 
range power law asymptotic behavior and diffusion processes can assume anomalous 
character whereas escape from a repeller may be algebraic rather than exponential. 
In quantum mechanical systems linked to a classical deterministic dynamics via 
periodic orbit trace formula, intermittency can be shown to lead to localisation and 
quasi-regular eigenfunctions 12 ) ~~ 14 ) . 

The Ruelle thermodynamical formalism 15 ) -17 ) relates dynamical averages to 
the spectral properties of the Perron-Frobenius operator or generalizations thereof, 
the Ruelle-Araki operators. For a 1-dimensional map a Perron-Frobenius operator 
relates the density at current time to the density at the previous time step: 

p n+1 (x) = C[ Pn ]{x) = J dy5(x - f{y)) Pn {y) = £ \f (y)^ p n (y)- (1-4) 

Here we shall show how the method of cycle expansions may be applied to the 
calculation of spectra of Perron-Frobenius operators in the presence of marginally 
stable orbits. We start in sect. 2.1 by recapitulating the theory for hyperbolic maps, 
and introduce dynamical zeta functions, spectral determinants, and their cycle ex- 
pansions. In sect. 2.2 we describe Perron-Frobenius operators in the presence of 
isolated marginally stable fixed points and introduce the induced map of an inter- 
mittent map. The core of the paper is sect. 3; we explain how marginally stable 
orbits are incorporated in cycle expansions. In order to gain some appreciation for 
the difficulties lying ahead, we start in sect. 3.1 by studying a piecewise linear map 
which follows the scaling (1-1) in the neighborhood of a fixed point, show that the 
dynamical zeta function has a branch cut, and extract from this cut the power law 
that characterizes the escape from the repeller. 

Such branch cuts are typical also for smooth intermittent maps with isolated 
marginally stable fixed points and cycles. In sect. 4 we discuss higher order curva- 
ture corrections required by their cycle expansions. The knowledge of the type of 
singularity one encounters for such systems enables us to construct the resummation 
method presented in sect. 4.1. Analytic continuations of dynamical zeta functions 
in terms of their integral representations are discussed in sect. 4.2. 

§2. Perron-Frobenius operators and cycle expansions 

Even though the operator formalism is not the focus of this paper, it is instructive 
to review it here as it sheds light on the cycle expansion approach pursued later. 
Briefly, when a dynamical system is uniformly hyperbolic and its symbolic dynamics 
is a subshift of finite type (a complete Smale horseshoe), spectral determinants of the 
associated evolution operators (like the Perron-Frobenius operator (1-4)) are entire 
and their spectra are discrete 18 ), very much like the bound spectra of quantum, 
electrodynamical and elastodynamic systems. However, if either the requirement of 
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hyperbolicity, or of the finite subshift, or both are not fulfilled, the spectral deter- 
minant is no longer an entire function, with delicate consequences for the spectra of 



We are far from having a global theory for non-hyperbolic dynamics; partial 
answers can, however, be given relating isolated features in the dynamics to spec- 
tral properties of the Perron-Frobenius operator. We will focus on a special case 
of non-hyperbolicity in what follows, that is, intermittent maps with complete sym- 
bolic dynamics and a single marginal fixed point. For intermittent dynamics where 
the symbolic dynamics is either not known, or cannot be finitely specified, a good 
practical method is the stability ordered cycle expansions 20 ) . Alternatively, one can 
make use of statistical methods developed by Dahlqvist 21 )' 22 ) in which spectral det- 
erminants are approximately described in terms of probability distributions for the 
return time to the non-hyperbolic region. 

Another example of non-hyperbolicity in the finite subshift setting is the Ulam 
map or logistic map x — > Axil — x), with strongly contracting non-hyperbolic critical 
point x c = 1/2. The Ulam map is conjugate to the tent map and this conjugacy 
fails only at the leftmost fixed point, the forward image of the critical point. It 
was shown in ref. 9) that for maps, such as the "skew" Ulam map which cannot be 
conjugated to a tent map, the spectrum can be obtained by dividing out from the 
spectral determinant the fixed point that the critical point iterates to. 

A similar technique needs to be employed when dealing with intermittency. As 
will be discussed later, marginally stable cycles are again factored out from the spec- 
tral determinant leading in this case, however, to a whole sub-determinant related 
to the dynamics close to the marginally stable point. The corresponding operator 
is (as can be shown by explicit conjugation) equivalent to free translation with a 
continuous spectrum, just as it is in quantum scattering where asymptotic states are 
also translation invariant. What remains of the factored spectral determinant can 
be interpreted as a spectral determinant of another dynamical system, the spectral 
determinant of a purely hyperbolic induced map. 

2.1. The hyperbolic case - a review 

In what follows we shall illustrate the key techniques by considering the dynamics 
of 1-dimensional hyperbolic maps x *— > f(x) on the unit interval M = [0, 1] with two 
piecewise-analytic monotone branches 



The two branches are assumed complete, /o(A'fo) = /l(-Mi) = M.. Such a map 
allows for escape if a < b, and the dynamics is bounded to the unit interval if 
a = b. Its symbolic dynamics is the full binary shift, with a 1-to-l mapping between 
admissible orbits and all possible infinite binary itineraries. A map is hyperbolic if 
the stability exponents of all periodic orbits are strictly positive 
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where the sum is over the cycle points Xi along the cycle p. In this case a well 
established theory exists 15 ^' 18 ), and it can be shown that the Perron- Frobenius oper- 
ators have discrete spectra {Ao, Ai, A2 . . .} when acting on a suitable space of real 
analytic functions. The method of cycle expansions 2 ^ ' 24 )> 17 ) is often an efficient 
method of computing the spectrum. It makes use of the fact that the trace of the 
Perron-Frobenius operator can be written in terms of the periodic orbits of the map, 



,1 W J 

tvC n = J2K= dxS(x-r(x)) = J2T T —rr < 2 " 3 ) 

JO ■ I 1 yl «l 



where the sum is taken over all periodic points of period n. The spectral determinant 
of £ is related to the traces by the cumulant expansion 17 ) , 

det (1 - zC) = 1 - ztuC - ^z 2 [trC 2 - (tr£) 2 ] - . . . 

=nn(i-pOT )=nff'w, <^> 

fc=0 p V I PI p/ fc=0 

where the product runs over all prime periodic orbits (repetitions excluded) of periods 
n p . For hyperbolic maps with complete binary symbolic dynamics, the spectral det- 
erminant is an entire function 15 )' 18 ) whose zeros {zq, z±, Z2, ■ ■ ■} are related to the 
eigenvalues of the Perron-Frobenius operator via A Q = z^. 1 . 

A quantity easily accessible to experiment and numerical calculations is the 
survival probability P n , that is, the fraction of phase space which remains within the 
phase space volume M (with M = [0, 1] for the map (2-1)) for at least n iterations 

Each surviving interval Mi contains a periodic point X{ £ Mi belonging to a prime 
cycle p or its rth repeat, of total period rn p = n. For hyperbolic maps their size can 
be asymptotically bounded from below and above in terms of the stabilities of that 
orbit, 

C» <J^!<C * (2-6) 
\A P \ \M\ \A P \ 

with Ci, C2 independent of p. The asymptotic behavior of P n is thus captured by the 
periodic orbit sum 

w 1 

f«~r„ = Y, w (2-7) 

For strictly hyperbolic maps this survival probability falls off exponentially with n 
in the large n limit. The periodic orbit sum (2-7) can be rewritten in terms of 1/Co 
defined in (2-4), 
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where the contour encircles the origin in the clockwise direction and lies inside the 
unit circle \z\ = 1. The dynamical zeta function l/£ can be written as a convergent 
product over periodic orbits in this region, integrals and sums can be interchanged, 
and the formula (2-7) is recovered after evaluating the integrals term by term by 
Cauchy integration. For hyperbolic maps, cycle expansions (discussed in more detail 
below) provide an analytic extension of the dynamical zeta function beyond the 
leading zero; we may then deform the original contour into a larger circle with 
radius R which encircles both poles and zeros of 1/Co- Residue calculus turns this 
into a sum over the zeros z a and poles zp of the dynamical zeta function, 

zeros poles 
\z a \<R 01 \z \<R P J ^R 

where the last term gives a contribution from a large circle 7^. We thus find ex- 
ponential decay of F n dominated in general by the leading zero zq of l/( (z). We 
expect zq = 1 for bounded dynamics and \zq\ > 1 for maps which allow for escape. 

The all important caveat here is whether one is able to construct an analytic 
continuation of l/( (z) for \z\ > \zo\. For systems with a complete symbolic dynamics 
with finite alphabet the method of cycle expansions 9 ) , with cycles and products of 
cycles (pseudo-cycles) ordered by their total period, has been demonstrated to work 
very well. For maps with a binary symbolic dynamics one has 

/ Tip \ 00 

1 /c w = n( 1 -rn) = 1 -E c « ( 2 - 10 ) 

1 n=l 



1 - Z 



1 1 

+ 



LIA)| |4| 



z 2 



1 



_ I |44| 



z 3 [■ • •] 



Here the c\ term, the fundamental term, includes all unbalanced, not shadowed 
cycles. For the higher terms, called curvature corrections, large cancellations between 
cycles and pseudo-cycles that shadow them take place 9 )' 17 ). 

For a piecewise linear map like the tent map (1-3) with \A p \ = 2 np for all periodic 
orbits of length n p , contributions from cycles and pseudo-cycles cancel exactly and 
the curvature contributions are identically zero. The dynamical zeta function thus 
has the form 

l/( Q (z) = l-z (2-11) 

and vanishes at zq = 1, as expected. In this case the complete spectrum of the Perron- 
Frobenius operator {Ao, Ai, A2, • • •} can easily be deduced from a cycle expansion of 
higher order dynamical zeta functions 1/Q k (z)\ one obtains X a = z^ 1 = l/2 2a , 
a = 0,1,2,.... 

2.2. The intermittent case and the induced map 

For intermittent maps of the form (IT) the trace and determinant formulas (2-3), 
(2-4) have no meaning as they stand; the periodic orbit stabilities can be arbitrarily 
close or equal to 1, and that renders the formulas singular. A spectral analysis of 
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the Perron-Frobenius operator is now much subtler than in the hyperbolic case. In 
this section we briefly sketch the main ideas behind the recent breakthroughs in 
the spectral analysis for intermittent maps 25 ) -28 ) which establish that the Perron- 
Frobenius operator has a continuous spectrum o~ c = [0, 1], and a pure point spectrum 
consisting of isolated eigenvalues of finite multiplicity. Due to the symmetry of the 
map in this case, zero is an eigenvalue of infinite multiplicity. We shall return to 
dynamical zeta functions and spectral determinants in sect. 3.2 and will show there 
that marginal orbits induce branch cuts in zeta functions which are related to the 
continuous part of the spectrum. 

The starting point for a rigorous analysis is the precise specification of the func- 
tion space that the Perron-Frobenius operator acts on. This requires a level of 
analysis that lies beyond the scope of this paper: in particular, one needs to prove 
that L maps this function space into itself. In what follows we will only give the key 
steps of the operator approach, which will be instructive later when discussing cycle 
expansions, without attempting to give details of the proofs. 

The basic idea is to separate the intermittent part of the dynamics from the 
chaotic dynamics. For clarity, we shall illustrate the method using the Farey map 
(1-2) as an example 27 ^' 28 \ Denote by Fi(x), i € {0, 1} the inverse of the branch /j 
with Fq(x) = x/(x + 1) and F±(x) = l/(x + 1). maps the unit interval A4 into 
Mi, the interval supporting fi. For the Farey map, the Perron-Frobenius operator 
(1-4) acts on the two branches as 

£ = Co + £i , 

£ »M«=(TT^'(^)' ^"(TW'G^t)' (2 ' 12) 

We focus on the intermittent behavior of the - branch by formally factoring it out 

(1 - zC) = (1 - t{z)) ■ (1 - z£o) ■ (2-13) 

Expanding 

(1 - t(z)) = (1 - z£o - *— - = 1 - zd—^-— 

1 — ZLq 1 — ZLq 

as a geometric series yields a formula for C: 

oo 

C(z) = * m ^- X ■ (2-14) 

m=l 

Loosely speaking, the partial Perron-Frobenius operator Cq in the decomposition 
(2-13) is responsible for the continuous part of the spectrum, a new feature specific 
to intermittent dynamics, while C adds isolated eigenvalues. 

The operator Cq is conjugated to the shift operator S[(p](x) = ip(l + x) under 
the transformation p{x) = C[(p\{x) = (p{fi{x)). The shift operator has a continu- 
ous spectrum a c = [0, 1] on a function space of analytic functions obtained by a 
generalized Laplace transform 

1 t'°° 

c^(x) = F[1>]{x) = — / dfi(s) e- s / x e- s i>{s) 
x Jo 
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acting on the space of square integrable functions ip £ L 2 (IR + ,/i) with dfi(s) = 
s/(e s — 1) ds 27 ) . Similar relations are expected to hold in neighborhoods of marginally 
stable fixed points for general intermittent maps 25 -*. 

It remains to understand the role of C. The meaning of the formal power series 
(2-14) is best grasped by setting z = 1: the operator 



£(1) = t=Y j C 1 C 



m— 1 




m=l 



is the Perron- Frobenius operator for the map 9 ** 



/(*) 



' fi(x) x £ M\ 

f2(x) = fi(fo(x)) x £ M 2 = 

' ; , (2-15) 

fm{x) = h{fo {x)) X G M m = (Qm, Qm-l] 



constructed from compositions of the initial 2-branch map / defined in (2T). The 
interval boundaries q n are successive preimages fo(q n ) = q n -i of the rightmost point 
of M.q, that is, q\ = 1/2 for the Farey map or q\ = a in case of the map (2-1). The 
map / thus has an infinity of branches, taking care of all sequences of consecutive 
iterations of the left branch of / in a single step. The length of a cycle of / thus 
equals the number of l's in the binary itinerary of the original map /. For the Farey 
map the map (2-15) is the Gauss map of number theory 10 ) 

/Gauss (x) = |^| x € (0, 1], 

where {• • •} denotes the fractional part. This map is related by conjugacy to the 
induced map of ref. 29), 30), the first return map on interval [1/2, 1]. The key issue is 
that the induced map is hyperbolic, so the operator C adds only isolated eigenvalues 
to the spectrum. 

Now, factoring out parts of the spectral determinant implicit in (2-13) does not 
necessarily make sense. If the spectral determinant is an entire function, as it is for 
hyperbolic maps, this merely introduces singularities in the remainder term and the 
spectrum of Cq may have little in common with the spectrum of the full Perron- 
Frobenius operator C Things are, however, different for intermittent maps. The 
determinant det (1 — C(z)) is analytic in C — [l,oo), with a branch cut along the 
real line [l,oo), and can be analytically continued through this branch cut, albeit 
onto different Riemann sheets. Hence a c = [0, 1] is part of the spectrum of the full 
operator C, while the rest of the spectrum of C is discrete and can be deduced from 
C. From (2-13) it follows now that z^ 1 is an eigenvalue of C if 1 is an eigenvalue of 
C(z). 

Let us finally return to the spectral determinant of C. The dynamics of / and 
/ are equivalent apart from singling out the fixed point at x = 0. As a consequence, 
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very little changes when writing the spectral determinant (2-4) for C(z) in terms of 
periodic orbits except for taking out the marginal fixed point, that is, 



det (1-4.) =nn(i-pOTj> p' i6 > 

where n p is again the length of the periodic orbit with respect to the original map 
/. That something more profound is happening becomes clear when we consider the 
first few terms in the cumulant expansion of the spectral determinant Eqn. (2-4), 

det (l - C{zj) = 1 - tr t(z) - . . . = 1 - f; — 1 - • • • • (2-17) 

m=l ' 

The infinite sum over orbits 10 m_1 reflects the fact that the induced map has in- 
finitely many fixed points. 



§3. Periodic orbit expansions for intermittent maps 



In order to get to know the kind of problems which arise when studying dynam- 
ical zeta functions in the presence of marginal stability we shall train on a carefully 
crafted piecewise linear map first, and return to the case of smooth intermittent 
maps in sect. 4. 

3.1. A piecewise linear toy map 

The tent map (1-3) is an idealized example of a hyperbolic map. To construct 
a piecewise linear analogue for intermittency, we start from the induced map (2T5), 
linear in each branch 31 ). 



Fig. 2. A piecewise linear intermittent map 
of type (2-1) for intervals chosen accord- 
ing to (3-5) with a = .5, b = .6, s = 1.0. 




We take the right branch in (2-1) expanding and linear: 

fl{x) = Y~ri^ x ~ 6 ) for x ^ M i = (P, 1] • 

The left branch /o of / is constructed in such a way that it will exhibit the intermit- 
tent behavior (IT) near the origin. Motivated by the form of the induced map (2-15), 
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we consider a monotonically decreasing sequence q m of preimages of the rightmost 
point of fo in Mq = [0, a], with q\ = a and q m — > as m — > oo. This sequence 
defines a partition of the left interval A4q into an infinite number of monotonically 
decreasing intervals Ai m , m > 2 with 

oo 

M m = (q m ,qm-i] and M = [j M m . (3-1) 

m=2 

Conversely, any given monotone sequence {q m } defines a piecewise linear con- 
tinuous map fo(x), with the slopes of the various linear segments given by 

fi( x ) = ih forxeJWi 

fo( x ) = forx€X 2 ( 3 . 2 ) 

foW = ifef for x € Mm > m ~ 3 



with |At m | = g m -i - q m for m > 1. 

We have seen in (2-15) and (2-17) that the family of periodic orbits with itinerary 
10 m_1 plays a key role for intermittent maps of the form (IT). An orbit 10 m_1 enters 
the intervals M.\ — > A4 m — > A4 m ~i, ■ ■ ■ — ► successively and the family comes 
closer and closer to the marginally stable fixed point at x = as m — > oo. The 
stability of a cycle 10 m_1 for m > 2 is proportional to the size of the interval closest 
to the marginal fixed point 

A 10 m-i = /o(x m )/o(x m _i) . . . fd(x 2 )f{{xi) = ^ ' ( 3 ' 3 ) 

with Xi € Mi- 

A piecewise linear map of the form (2T) is fixed by the sequence {q m }- By 
choosing q m = 2~ m , for example, we recover the uniformly hyperbolic Bernoulli shift 
map. An algebraically decaying sequence {q m } 

' \M m \~ -7^777, (3-4) 



Mm i/i l'"ra| 1-1-1 /e ' 

yields an intermittent map of the form (3T), with the asymptotic behavior controlled 
by the intermittency exponent s in (IT). The stability eigenvalues of periodic orbit 
families approaching the marginally stable fixed point, such as the 10 m_1 family in 
(3-3), grow in turn only algebraically (not exponentially) with the cycle length. 

It may now seem natural to construct an intermittent toy map in terms of a 
partition \A4 m \ = l/m 1+1 / s , that is, a partition which follows (3-4) exactly. Such a 
choice leads to a dynamical zeta function which can be written in terms of Jonquiere 
functions (or polylogarithms) 32 )' 33 ) . We shall, however, use the freedom in choosing 
the partition to make life as simple as possible later on, without loosing the key 
intermittency features. Inspired by a bit of reverse engineering we fix the intermittent 
map intervals M. m by 

„r(m + £-l/s-l) „ , , 
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where -T(x) is the Gamma function, £ = [1/s] denotes the integer part of 1/s, and 
the normalization constant C is fixed by the full partition condition Ylm=2 \-M m \ = 
q\ = \Mo\, that is, 

rim -lis)' 



C= \M \ 



E 

.m=l+l 



r(m + 1) 



(3-6) 



One easily verifies that the intervals decay asymptotically like required 
by the condition (3-4) using Stirling's formula for the Gamma function. 

The point spectrum of / is given by the zeros of the spectral determinant (2-14) 
or the associated dynamical zeta functions. In order to find analytic continuations 
of the dynamical zeta function l/C^j we start by expanding the product as we did 
in (2-10). 

Due to the absence of the fixed point 0, there are no pseudo-cycles that shadow 
cycles of the form 10 m_1 , so the fundamental cycles are now the entire infinite family 
of cycles with itineraries of the form 10 m_1 accumulating towards the marginal fixed 
point 0. As the fixed point does not participate in the dynamics, we have to switch 
from the 2-letter {0, 1} alphabet to the infinite alphabet m = 10 m_1 which labels the 
branches of the induced map (2-15). The admissible itineraries are now generated 
by all walks on the infinite Markov diagram 9 ) fig. 3. 



Fig. 3. Markov graph for the infinite alpha- 
bet {m = 10 m_1 ;0, m > 1} derived 
from full binary symbolic dynamics by 
separating out the fixed point. 




The piecewise linear form of the map which maps intervals A4 m exactly onto 
M. m -\ simplifies matters considerably: when expanding the zeta functions in the 
spectral determinant (2-16) all cycles which traverse the right branch at least twice 
are cancelled exactly by pseudo-cycles. The cycle expanded dynamical zeta functions 
depend thus only on the fundamental cycles and have the form 



(3-7) 



^ |/l 10 m-l |/l^ Qm _ 1 

To keep the discussion simple, we will restrict our considerations to 1/Co- We obtain 



, . , , , , x „ 1 — 6 rim + 

l/C (z) = l-(l-6)z-C— - 

m=2 



r{m + 



(3-8) 



The cycle expansion does not, however, lead to an analytic continuation here, despite 
the exact cancellation between cycles and pseudo-cycles, in contrast to the dynam- 
ical zeta function for the tent map (1-3). Due to the slow, algebraic decay of cycle 
weights, the sum (3-8) is divergent for \z\ > 1. We shall now show that this behavior 
is due to a branch cut in l/( k (z) along [1, oo). 
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3.2. Branch cuts in dynamical zeta functions due to intermittency 

In order to analytically continue the sum (3-8) beyond \z\ > 1 we make use of 
the following binomial identities valid for a = l/s > 0: 

a non-integer (1 - z) a = ^ - ^-^— z n (3-9) 

^-^ 1 \—ol)1 (n + 1) 

(l-zriog(i-z) = £:=i(-i) n c„z n 

"integer +(_iy*+i n ,!v 00 -n 



v 7 fc=0 



For simplicity sake, we restrict the intermittency parameter to 1 < l/s < 2, that is 
the case [l/s] = £ = 1. All that follows can easily be generalized to arbitrary s > 
using (3-9). What once was a mystifying choice of interval widths (3-5) now pays 
off: the infinite sum (3-8) can be evaluated analytically, 

^ r(m-l/s) / r(-I)[(l-z)V-_i + I z ] for 1< l/s < 2; 

^ r(m + 1) " \ (1 - z) log(l - *) + z for s = 1 . 

The normalization constant C in (3-6) can be evaluated, and one obtains for 1 < 
l/s < 2: 

1/CoW = 1 - (1 - &)z - 1/g a _ 11 1 I^ ((1 - ^) 1/s - 1 + ^) (340) 



and for s = 1: 



I/CO*) = 1 - (1 - 6)z - a|— ^ ((1 - z) log(l - z) + z) . (3-11) 



For general s > one obtains 

« L - ft I 

ft 

for non-integer s with £ = [l/s] and 



1/CoW = 1 - (1 - b)z ~ r^n^^T ((1 " *) 1/S ~ 9s{z)) (3-12) 



1/C (z) = 1 - (1 - &)* - ^1) r^^T ((! " ^ lo §(l - *) ~ 9s(z)) (3-13) 

for l/s = £ integer; here, g s (z) are polynomials of order £ = [l/s] which can be 
deduced from (3-9). We see that the dynamical zeta function has a branch cut 
starting at z = 1 and running along the positive real axis. We find algebraic branch 
cuts for non integer intermittency exponents l/s and logarithmic branch cuts for 
l/s, integer. We will see in sect. 4 that branch cuts of that form are generic for 
1-dimensional intermittent maps. In the next section, we will use these analytic 
expressions to calculate the escape from our toy map for a < b in (2-1). 
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3.3. Escape rate 

The fraction of survivors r n defined in (2-5) have for hyperbolic maps been linked 
to 1/Co w hh the help of bounds given in (2-6). These bounds are no longer valid for 
intermittent maps. It turns out, however, that bounding survival probabilities strip 
by strip is a far too strict requirement for establishing the relation (2-7). In fact, 
a somewhat weaker bound can be established, linking the average size of intervals 
along a periodic orbit to the stability of the periodic orbit for all but the interval 
A4o«. We can write this bound by averaging over each prime cycle p separately, 



Cy 



1 



A 



< 



1 



n,. 



£ 

p 



\M^\ 
\M\ 



1 



I A 



(3-14) 



for some positive constants Cy C2 independent of p. A proof, which relies on the 
hyperbolicity of the induced map, can be found in ref. 34). Summing over all periodic 
orbits leads then again to (2-7) which is linked to 1/Cn via the integral representation 
(2-8). 



Fig. 4. The survival probability _T„ calcu- 
lated by contour integration; integrating 
(2-8) inside the domain of convergence 
\z\ < 1 (shaded area) of l/£ (z) in peri- 
odic orbit representation yields (2-7). A 
deformation of the contour j~ (dashed 
line) to a larger circle 7^ (without cross- 
ing the branch cut) gives contributions 
from the poles and zeros (x) of l/£ (z) 
between the two circles, as well as con- 
tributions from the integration along the 
branch cut 7 

C U t • 

The path deformation to \z\ > 1 that led to (2-9) requires more care here, as it 
must not cross the branch cut along [1, 00]. When expanding the contour to large 
\z\ values, we have to deform it in such a way that it sandwiches the branch cut in 
anti-clockwise direction, see fig. 4. Denoting the detour around the cut as Jcut, we 
may write symbolically 

zeros poles 

£-£ + 



Im z 




V 


® \ 


--J--J-J--J4'- 


1 < - — 


L a 


i_' z=l lit- . 1 




: ® J 







7h 



7cut 



where the sums include only the zeros and the poles in the area enclosed by the 
contours. 

Let us now go back to our intermittent toy map. The asymptotics of the survival 
probability of the map is governed by the behavior of the integrand 4- log Cq 1 in (2-8) 
at the branch point z = 1. We restrict ourselves again to the case 1 < 1/s < 2 first 
and write the dynamical zeta function (3T0) in the form 

1/C„(*) =ao + ai(l-z) + b (l - zfl* = G(l - z) 
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with 



Setting u = 1 



ao = z , 

1 — a 

z, we need to evaluate 



1- 1/sl 



1 

2vri 



[I 



"I cut 



du 



log G(u)du 



(3-15) 



where j cu t goes around the cut (that is, the negative u axis). Expanding the inte- 
grand -4- log G(u) = G'(u)/G(u) in powers of u and u l / s at u = 0, one obtains 



r/ 



logG(u) 



ai 1 6n 

— + — -'. 



+ 0{u) . 



du " ' ao s ao 
The integral along the cut may be evaluated using the general formula 

T(n-a-l) 1 



1 

2ri 



u a {l-u)~ n du 



Jcut 



r(n)r(—a) n 



a+l 



(l + 0(l/n)) 



(3-16) 



(3-17) 



which can be obtained by deforming the contour back to a loop around the point 
u = 1, now in positive (anti-clockwise) direction. The contour integral then picks up 
the (n — l)-st term in the Taylor expansion of the function u a at u = 1, see (3-9). 
Inserting (3T6) into (3-15) and using (3-17) we get the asymptotic result 



i 



i 



i 



i 



i 



M 

a s r(l - l/s) n l / s s - 1 b - a T(l - 1/s) n 1 / 8 



(3-18) 



We see that, asymptotically, the escape from an intermittent repeller is described by 
power law decay rather than the exponential decay we are familiar with for hyperbolic 
maps. The zeros and poles of 1/Co & ve exponentially decreasing contributions, which 
may indeed be dominant over large time interval before the power law decay sets in 
for n — > oo. A numerical simulation of the power-law escape from an intermittent 
repeller is shown in fig. 5. 




Fig. 5. The escape from an intermittent re- 
peller typically starts out exponentially, 
controlled by the zeros close to the cut 
but beyond the branch point z = 1, as 
in fig. 4, followed by a power law decay 
controlled by the type of the branch cut. 



10 




200 400 



600 800 



For general, non-integer l/s > 0, we write 

l/C Q {z) = A{u) + {u) 1 / s B{u)^G(u) 
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with u = 1 — z and A(u), B(u) are functions analytic in a disc of radius 1 around 
u = 0. The leading terms in the Taylor series expansions of A (u) and B{u) are 



b — a a 1 

ao = , »o 



g s (l) 1-a' 

see (3T2). Expanding -^-logCr('u) around u = 0, one again obtains leading order 
contributions according to (3T6) and the general result follows immediately using 
(3-17) 

a 1-6 11 
n ~ sg s (l)b-ar(l-l/s)n^- (3 ' 19) 

Applying the same arguments for integer intermittency exponents 1/s = £, one 
obtains 

r n ~ ( _1 ) <+1 — TTTI— "~f- P-20) 
sge(l) o — an 1 

So far, we have considered the survival probability for a repeller, that is we 
assumed a < b. The formulas (3-19) and (3-20) do obviously not apply for the 
case a = b of a bounded map. The coefficient ao = (b — a)/(l — a) in the series 
representation of G(u) is zero and we need to take into account next leading terms 
in the expansion of the logarithmic derivative of G(u) in (3-16). One obtains 



du 



**G(«) = ( I l + o^-V.) ,>i 



where we assume 1/s non-integer for convenience. The survival probability is thus 
in leading order 

\ + 0{n l - l ' s ) 8 <1 
l/s + Oin 1 / 3 ' 1 ) 8 >1 ' 



r n ~ 



For s < 1, this is what we expect. There is no escape, so the survival probability 
is equal to 1, which we get as an asymptotic result here. The result for s > 1 is 
somewhat more worrying. It says that F n defined as sum over the instabilities of the 
periodic orbits (2-7) does not tend to unity for large n. However, the case s > 1 is in 
many senses anomalous. For instance, the invariant density cannot be normalized 35 ) . 
It is therefore not reasonable to expect that periodic orbit theories will work without 
complications. 



§4. Smooth intermittent maps 



Now we turn to the problem of constructing cycle expansions for general smooth 
1-dimensional maps with a single isolated marginal fixed point. To keep the notation 
simple, we will consider two-branch maps with a complete binary symbolic dynamics 
as in (2-1). The necessary modifications for multi-branch maps with an isolated 
marginal fixed point are straightforward. We again assume that the behavior near 
the fixed point is given by (IT). This implies that the stability of a family of cycles 
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approaching the marginally stable orbit, for example the family 10 m , will increase 
only algebraically. Again we find for large m 

1 1 

AiQm m} +l l s ' 

where s denotes the intermittency exponent. 

Zeta functions or spectral determinants in terms of products over cycles are 
formally identically to those of hyperbolic systems except for the omission of the 
marginal orbit 0, see (2-16). When considering cycle expansion as in sect. 2.1, pe- 
riodic orbit contributions from, for example, cycles with itinerary 10 m are unbal- 
anced and we arrive at a cycle expansion in terms of infinitely many fundamental 
terms 13 ** as for our toy map. This corresponds to moving from our binary symbolic 
dynamics to an infinite symbolic dynamics by making the identification 

lO"^ 1 -> m; 10 fc_1 10 m - 1 -> km; 10 , - 1 10 fc-1 10 m_1 -» Ikm; 

table I. The topological length of the orbit is thus no longer determined by the 
iterations of the two-branch map /, but by the number of times the cycle reaches 
the right branch, that is, by the number of symbols 1 in the itinerary, or by the 
number of letters in the symbolic dynamics of the induced map / of /. 



oo - alphabet 






binary alphabet 








m = 1 


m = 2 


m = 3 


m = 4 


m = 5 


1-cycle 


m 


1 


10 


100 


1000 


10000 


2- cycle 


km 














lm 


11 


110 


1100 


11000 


110000 




2m 


101 


0101 


10100 


101000 


1010000 




3m 


1001 


10010 


100100 


1001000 


10010000 




4m 


10001 


100010 


1000100 


10001000 


100010000 


3- cycle 


Ikm 














11m 


111 


1110 


11100 


111000 


1110000 




12m 


1101 


11010 


110100 


1101000 


11010000 




13m 


11001 


110010 


1100100 


11001000 


110010000 




21m 


1011 


10110 


101100 


1011000 


10110000 




22m 


10101 


101010 


1010100 


10101000 


101010000 




23m 


101001 


1010010 


10100100 


101001000 


1010010000 




31m 


10011 


100110 


1001100 


10011000 


100110000 




32m 


100101 


1001010 


10010100 


100101000 


1001010000 




33m 


1001001 


10010010 


100100100 


1001001000 


10010010000 



Table I. Infinite alphabet versus the original binary alphabet for the shortest periodic orbit families. 

For generic intermittent maps, curvature contributions in the cycle expanded 
zeta function will not vanish exactly as they do for our piecewise linear model. 
The most natural way to organize the cycle expansion in this case is to collect 
cycles and pseudo-cycles of the same topological length with respect to the infinite 
alphabet of the induced map. Denoting cycle weights in the new alphabet as tfc m ... = 



Cycle expansions for intermittent maps 



17 



t 10 fc-iio™-i...i one obtains 13 ^ 



oo 



co _i =ii( i - t f) =i -E c ™ ^ 

p^O n=l 



oo oo oo 



— 1 ^ ] t m ^ ^ ^ ] 2 (^fe™ t^tm) 
m=l k=l m=l 

~ Y Y Y [ 3 tikm ~ 2 tiktm + Q titktm ) ~ Y Y Y Y ■ ■ ■ ■ 

1=1 k=lm=l ^ ' j=l 1=1 Hm=l 

with t p = z np /\A p \ and n p is the length of the binary itinerary. The first sum is the 
fundamental term, already noted in the cycle expansion (3-7) for the toy model. The 
curvature terms c n for n > 2 are now n-fold infinite sums where the prefactors take 
care of double counting of prime cycles. 

We consider the fundamental term first. For generic intermittent maps, we can 
not expect to obtain an analytic expression for infinite sums of the form 

oo 

f{z) = Y J Kz n (4-2) 

n=0 

with algebraically decreasing coefficients h n ~ n^^/ s+1 \ 

To evaluate the sum, we face the same problem as for our toy map: the power 
series diverges for \z\ > 1, that is, exactly in the 'interesting' region where poles, zeros 
or branch cuts of the zeta function are to be expected. By carefully subtracting the 
asymptotic behavior with the help of (3-9), one can in general construct an analytic 
continuation of f(z) around z = 1 of the form 36 ) 



f(z)~ A{z)+B{z){l- z) 1 ' 3 1/agN (4-3) 
f(z) ~ A(z) + B(z)(l - z) 1 / 3 ln(l -z) 1/s € N , 

where A(z) and B(z) are functions analytic in a disc around z = 1. From here we 
can derive the asymptotic behavior which coincides with the asymptotics obtained 
in (3T9) and (3-20), that is, the survival probability is in general given as 

r n ~ ~rr f° r n 00 (4-4) 

n l l s 

for 1-dimensional maps of the form (1-1). We have to work a bit harder if we want 
more detailed information like exponential precursors given by zeros or poles of the 
dynamical zeta function or higher order corrections. This information is buried in 
the functions A(z) and B(z) or more generally in the analytically continued zeta 
functions including curvature contributions. To get such an analytic continuation, 
one may follow two different strategies which we will sketch next. 
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4.1. Resummation 

One way to get information about the zeta function near the branch cut is to 
derive the leading coefficients in the Taylor series of the functions A(z) and B(z) 
in (4-3) at z = 1. This can be done in principle, if the coefficients h n in sums like 
(4-2) are known (as for our toy model). One then considers a resummation of the 
form ^ ~ ^ 

OO OO oo 

£ h S 2? = £ aj(l - zY + (1 - zf' s £ 6,(1 - zy, (4-5) 

3=0 3=0 3=0 

and the coefficients a,j and bj are obtained in terms of the hj's by expanding (1 — z) J 
and (1 — z)i +1 / s around z = using (3-9), and equating the coefficients. 

In practical calculations one often has only a finite number of coefficients hj, < 
i < un, which may have been obtained by finding periodic orbits and their stabilities 
numerically. One can still design a resummation scheme for the computation of the 
coefficients a,j and bj in (4-5). We replace the infinite sums in (4-5) by finite sums of 
increasing degrees n a and and require that 

n a n b n N 

^ai(l-zj + (\-z?l a ^bi{\-zj = Y / hi zi + °( znN+1 ) ■ ( 4 '6) 

i=0 i=0 i=0 

One proceeds again by expanding the right hand side around z = 0, skipping all 
powers z riN+1 and higher, and then equating coefficients. It is natural to require 
that |ra& + 1/s — n a \ < 1, so that the maximal powers of the two sums in (4-6) are 
adjacent. If one chooses n a + + 2 = + 1, then, for each cutoff length njy, 
the integers n a and n& are uniquely determined from a linear system of equations. 
The prize we pay is that the so obtained coefficients depend on the cutoff njy. One 
can now study convergence of the coefficients a,j, and bj, with respect to increasing 
values of un, or various quantities derived from a,j, and bj. Note that the leading 
coefficients ao and bo determine the prefactors as in (3-18). The resummed expression 
can also be used to compute zeros, inside or outside the radius of convergence of the 
cycle expansion ^ hjz 3 . 

The scheme outlined in this section tacitly assumes that a representation of the 
form (4-3) holds in a disc of radius 1 around z = 1. Convergence is improved further 
if additional information about the asymptotics of sums like (4-2) is used in the 
ansatz (4-5). 

4.2. Analytical continuation by integral transformations 

We will now introduce a method which provides an analytic continuation of 
sums of the form (4-2) without relying on the ansatz (4-5) 12 )' 13 ). The main idea is to 
rewrite the sum (4-2) as a sum over integrals with the help of Poisson summation and 
find an analytic continuation of each integral by contour deformation. In order to do 
so, we need to know the n dependence of the coefficients h n = h(n) explicitly for all 
n. If the coefficients are not known analytically, one may proceed by approximating 
the large n behavior in the form 

h(n) = n- l l s -\C 1 + C 2 n- 1 + ...), n^O, 
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and determine the constants Q numerically from periodic orbit data. By using the 
Poisson resummation identity 

oo oo 

y~] 5{x — n) = ^2 exp(27rimx) , (4-7) 

n=— oo m=—co 

we may write the sum (4-2) as 

-1 00 /»oo 

/(*) = -/i(0)+ V / dxe 2nimx h(x)z x . (4-8) 

The continuous variable x corresponds to the discrete summation index n and it is 
convenient to write z = rexp(ia) from now on. The integrals are of course still not 
convergent for r > 0. An analytic continuation can be found by considering the con- 
tour integral, where the contour goes out along the real axis, makes a quarter circle 
to either the positive or negative imaginary axis and goes back to zero. By letting 
the radius of the circle go to infinity, we essentially rotate the line of integration form 
the real onto the imaginary axis. For the m = term in (4-8), we transform x — > ix 
and the integral takes on the form 

/•oo roc 

/ dxh{x)r x e ixu = i dx h(ix)r ix e - xa . 
Jo Jo 

The integrand is now exponentially decreasing for all r > and a ^ or 2tt. The 
last condition reminds us again of the existence of a branch cut at Re z > 1. By the 
same technique, we find the analytic continuation for all the other integrals in (4-8). 
The real axis is then rotated according to x — > ±ix where ± sign refers to the sign 
of m, 

foo poo 

/ dxe ±2 ™\ m \ x h{x)r x e ixa = ±i / dx h(±ix) r ±^ e ^( 2 -H±-). 
Jo Jo 

Changing summation and integration, we can carry out the sum over \m\ explicitly 
and one finally obtains the standard formula 

f(z) = lh(0) +i f dx h(ix) r lx e~ xa (4-9) 
2 Jo 

poo —2nx 

+ i J dx - _ e _ 2wx [h(ix)r ix e- xa - h{-ix)r~ ix e xrj \ . 



The transformation from the original sum to the two integrals in (4-9) is exact for 
r < 1, and provides an analytic continuation for r > 0. The expression (4-9) is 
especially useful for an efficient numerical calculations of the dynamical zeta func- 
tion for \z\ > 1, which is essential when searching for zeros and poles of the zeta 
function. 
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4.3. Curvature contributions 

So far, we have discussed the fundamental term Ylri=i^n in (4-1) and showed 
ways to extract the information from such power series with algebraically decreas- 
ing coefficients. The fundamental term determines the main structure of the zeta 
function in terms of the leading order branch cut. Corrections to both the zeros and 
poles of the dynamical zeta function as well as the leading and sub-leading order 
terms in expansions like (4-3) are contained in the curvature terms in (4-1). The 
first curvature correction 

k=l m=l 

also has algebraically decaying coefficients which diverge for \z\ > 1. The analytically 
continued curvature terms have again branch cuts along the positive real z - axis. 
Our ability to calculate the higher order curvature terms depends on how much we 
know about the cycle weights tk m - The cycle expansion itself suggests that the terms 
tkm decrease asymptotically like 

~ < 4 ' 10 > 

for 2-cycles of the induced map and in general for n-cycles like 

1 

^mim2...m„ ~ 



m\vri2 • • • m n ) 1+1 / s 



If we know the cycle weights i mim2 ... mn analytically, we may proceed like in sect. 4.2 
by transforming the multiple sums into multiple integrals and by rotating the axis' 
of integration. 



§5. Summary and conclusions 



Marginally stable orbits affect the analytic structure of dynamical zeta functions 
and the rules for constructing cycle expansions: While the marginal orbits have to 
be omitted, the cycle expansions need to include families of infinitely many longer 
and longer unstable orbits which accumulate towards the marginally stable cycles. 
Correlations for such non-hyperbolic systems decay algebraically with the decay 
rates controlled by the branch cuts of dynamical zeta functions. Compared to pure 
hyperbolic systems, the physical consequences are drastic: exponential decays are 
replaced by slow power law decays, and transport properties, such as the diffusion, 
may become anomalous 40 )' 41 ). 
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